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ABSTRACT 

In this paper, we develop the method of analyzing the velocity field of cosmic 
matter with a multiresolution decomposition. This is necessary in calculating the 
redshift distortion of power spectrum in the discrete wavelet transform (DWT) 
representation. We show that, in the DWT analysis, the velocity field can be described 
by discrete variables, which are given by assignment of the number density and velocity 
into the DWT modes. These DWT variables are complete and not redundant. In 
this scheme, the peculiar velocity and pairwise velocity of galaxies or particles are 
given by field variables. As a consequence, the velocity dispersion (VD) and pairwise 
velocity dispersion (PVD) are no longer measured by number-counting or pair-counting 
statistic, but with the ensemble of the field variables, and therefore, they are free 
from the bias due to the number-counting and pair-counting. We analyzed the VD 
and PVD of the velocity fields given by the N-body simulation for models of the 
SCDM, rCDM and ACDM. The spectrum (scale-dependence) of the VD and PVD 
show that the length scale of the two-point correlation of the velocity field is as large 
as few tens h~^ Mpc. Although the VD and PVD show similar behavior in some 
aspects, they are substantially different from each other. The VD-to-PVD ratio shows 
the difference between the scale-dependencies of the VD and PVD. More prominent 
difference between the VD and PVD is shown by probability distribution function. 
The one-point distribution of peculiar velocity is approximately exponential, while the 
pairwise velocity's is lognormal, i.e. of long tail. This difference indicates that the 
cosmic velocity field is typically intermittent. 

Subject headings: cosmology: theory - large-scale structure of universe 



1. Introduction 

The radial (redshift) distance of galaxies is distorted by their peculiar velocity. Observation 
of redshift survey provides only a distorted picture of the galaxy spatial distribution in the radial 
direction. The galaxy power spectrum P (k) in the Fourier representation is measured from the 
redshift-distorted galaxy distribution. It is the power spectrum in redshift space. To recover the 
power spectrum in real space, P^{k), one needs to map the distributions of galaxies from the 
redshift space to the real space. Semi-phenomenological theory of redshift distortion show that 
the mapping between redshift and real spaces gives 

P^{k) = G{k)P^{k), (1) 

where redshift distortion factor G{k) depends on two parameters: redshift distortion parameter (5, 
and 1-D peculiar velocity dispersion (VD) a"" . This leads to two effects: 1. the enhancing of power 
on large scales due to the linear effect of redshift distortion; 2. the suppressing of power on small 
scales due to the random motions of galaxies inside virialized groups and clusters of galaxies. 

In deriving the relation (1), the following assumptions are employed: 1. the effect of the 
coupling between the peculiar velocity and the density perturbations is linear (Kaiser 1987); 2. 
the effect of the random motions of galaxies is independent of the density perturbations of cosmic 
mass field, 3. the probability distribution function (PDF) of random peculiar velocity or the 
pairwise velocity is either gaussian or exponential. Therefore, besides the linear effect, or point 1, 
eq.(l) does not require the information of the spatial distribution of peculiar velocity of galaxies 
(or dark matter). The parameter a^ actually is measured from the statistics of number-counting 
of galaxies or pair-counting of galaxies. That is, the spatial distribution of the peculiar velocity of 
galaxies or cosmic matter, v(x), is not treated as a random field. 

Recently, we developed the method of measuring the galaxy power spectrum with a space-scale 
(multiresolution) decomposition, i.e. measuring power spectrum in the representation of discrete 
wavelet transform (DWT) (Fang & Feng 2000, Yang et al. 2001). Unhke eq.(l), the redshift 
distortion of the DWT power spectrum must be described by the statistics of the velocity field 



v(x), as the number-counting or pair-counting statistics are not enough. For instance, the redshift 
distortion of the power of a DWT mode ip{x) is dependent on the projection of the velocity 
field on this mode, i.e. / v(x)^(x)(ix, which cannot be simply measured by number-counting or 
pair-counting of galaxies. 

The purpose of this paper is to develop the method of analyzing the velocity field of galaxies 
of dark matter particles as a random field. That is, we introduce the variables of the velocity 
field by an orthogonal DWT decomposition of the field, and calculate all statistics of velocity 
field by average over the ensemble of the field variables. With this field description, we study 
the problems referring eq.(l): 1.) the scale dependence of velocity dispersion (VD) and pairwise 
velocity dispersion (PVD); 2.) the PDF of peculiar velocity and pairwise velocity; 3.) the local 
correlation between the velocity field and density field of cosmic matter. 

Other motivation of studying velocity field is from the intermittency of cosmic density field. 
The analysis on high resolution data of QSO's Lya forests has revealed that the cosmic mass 
field is significantly intermittent (Jamkhedkar, Zhan &: Fang 2000; Feng, Pando & Fang 2001; 
Zhan, Jamkhedkar & Fang 2001). That is, the PDF of local density fiuctuations on small scales 
is neither gaussian nor exponential, but long-tailed. Taking into account the coupling between 
the local density fiuctuations and peculiar velocity, we can expect that the PDF of the velocity 
field should be neither gaussian nor exponential, but long-tailed on small scales. This study also 
requires a decomposition by orthogonal basis, as a superposition of independent variables may 
erase the long-tail feature due to the central limit theorem. 

The paper will be organized as follows. §2 introduces the multiresolution analysis of the 
velocity field. It will focus on the DWT description of the velocity dispersion (VD) and pairwise 
velocity dispersion (PVD). §3 develops the method of analyzing the scale-dependence of the VD 
and PVD, and demonstrate it with an uniformly random sample. §4 studies the VD and PVD of 
N-body simulated samples, including the scale dependence of the VD and PVD, the VD-to-PVD 
ratio, and the PDFs of the peculiar velocity and pairwise velocity. In §5, we show the correlation 
between the VD, PVD and the local density and local density fluctuations. Finally, the conclusions 



and discussions will be in §6. 

2. Velocity fields in the DWT representation 
2.1. Peculiar velocity dispersion and relative velocity dispersion 

For a given galaxy catalog, one can write down a number density distribution of galaxies as 

n(x) = ^it;„(5^(x-x„), (2) 

ra=l 

where Ng is the total number of galaxies in the catalog, x„ is the position of the nth galaxy, Wn 
is its weight, and 5^ the 3-D Dirac 5 function. n(x) is a sampling of the cosmic mass density 
distribution. 

The velocity field, v(x), describes the deviation of particle motion from Hubble flow. The 
peculiar velocity of galaxies, v(x„) and n = l...Ng, is a sampling of the velocity field. Thus the 
continuous velocity field v(x) can be constructed by the convolution of the particle velocity with 
an assignment function W{rj), 

^ / N ^ Y.n=i ^(|xn - x|/r)u;„v(x„) 
E^^ii^lxn-xl/rK 

where r is a filter scale. One of the simplest example of assignment functions is the top-hat 

window function, i.e. W{ri) = 1 for rj < 1/2 and for rj > 1/2. Thus, Vr(x) can be considered as a 

sampling of the cosmic velocity field v(x) smoothed on the scale r. 

Using Vr(x), the velocity dispersion (VD) is defined as the variance of the one-point 
distribution of the random variable Vj.(x), i.e. 

K(x)]2 = (v2(x)). (4) 

where (...) means the ensemble average. 

Using eq.(3), one can define the relative velocity by 

Av,(x) = v,/2[x + (r/2)] - v,/2[x - (r/2)], (5) 
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where |r| = r. Obviously, Avr(x) measures the difference of the velocity field at x + r/2 and 
X — r/2. If r is small, and there is only one galaxy x^ in the sphere r/2 at x + (r/2), and one 
galaxy x^ in the sphere r/2 at x — (r/2), we have Avr(x) = v(xn) — v(xm). Therefore, the 
measure eq.(5) contains all information of the relative velocity of individual galaxy pairs. 

Yet, the relative velocity defined by eq.(5) is not completely the same as pairwise velocity of 
galaxies. Usually, the pairwise velocity is defined as t;i2 — vu, where ^12 is the relative velocity 
of galaxy pair 1 and 2 along line-of-sight, and IJ12 is the mean relative velocity, or the in- fall 
motion, which generally is modeled by a similarity solution (Davis & Peebles 1983). Therefore, 
the measure of pairwise velocity is model-dependent. 

The reason of using the relative velocity eq.(5) mainly is that the variables Avr(x) are able 
to provide an information-lossless and model-independent decomposition of the velocity field 
(§2.2). Moreover, the in-fall term of the pairwise velocity measure is negligible for small scales. 
Therefore, at least on small scales, the definition eq.(5) is a reasonable generalization of the 
pairwise velocity of galaxies on small scales. We may still call Avj.(x) the variables of pairwise 
velocity decomposition. 

Thus, the pairwise velocity dispersion (PVD) is defined as the variance of the one-point 
distribution of the random variable Avr(x), i.e. 

[ar(x)]2 = ([Av,(x)]2). (6) 

2.2. Need for scale-space decomposition 

A common feature of v,.(x) [eq.(3)] and Av,.(x) [eq.(5)] is that each quantity depends on two 
variables: the position x and the scale r. That is, the VD and PVD measure the random velocity 
field under a space-scale decomposition. This point can be seen more clearly by rewriting eqs.(3) 
and (5) using eq.(2). We have 

^ /n(xOv(xOT^(|x--x|/r)dx- 
"■^ ' /n(x')VF(|x'-x|/r)dx' ' ^ ' 



and 



Avrfx) 



/n(x')v(x')[/r(x'-x)dx' 



/n(x')VF(|x'-x|/r)dx' ' 
where function U{x' — x) is given by the subtractions of window functions, i.e. 

Ur{^' - x) = VF[|x' - X + r/2|/(r/2)] - VF[|x' - x - r/2|/(r/2)], 



(8) 



(9) 



Ur (X — X) 



or 

1 |x'-x-r/2| < r/2 

-1 |x' -x + r/2| < r/2 
otherwise. 
Function [/^(x — x') is similar to the Haar wavelet, of which in 1-D is 



V'(x) 



1 < X < 1/2 
-1 l/2<x<l 
otherwise. 



(10) 



(11) 



Therefore, the pairwise velocity eq.(8) actually is a decomposition of the velocity field with a Haar 
wavelet-like function. 

However, the functions Ur{x' — x) with respect to x and r don't exhibit completeness and 
orthogonality. Therefore, the decomposition with [/^(x' — x) may lead to loss of information if 
C/'s are incomplete, or cause false correlations if they are redundant. To have a proper measure 
of the VD and PVD, we call on discrete wavelet transform (DWT), which provides a complete 
and unredundant space-scale decomposition. The first approach of describing cosmic velocity field 
by wavelet is given by Rauzy, Lachieze-Rey & Henriksen (1993). However, they use continuous 
wavelets, which give rise to a redundant decomposition (Farge 1992). 



2.3. The DWT variables of velocity field 

For the details of the mathematical properties of the DWT refers to Mallat ( 1989a, b); Meyer 
(1992); Daubechies, (1992), and for physical applications, refer to Fang & Thews (1998). The 



Haar wavelet [eq.(12)] provides a clear picture of the DWT decomposition, and it is also easy for 
numerical work. However, the Haar wavelet is discontinues, and therefore, it is not well localized 
in scale space. To our work, the most important properties of the basis for the scale-space 
decomposition are 1.) orthogonality, 2.) completeness, and 3.) locality in both scale and physical 
spaces. Therefore, wavelets with compactly supported basis are suitable to the velocity field 
analysis. Among the compactly supported orthogonal basis, the Daubechies 4 (D4) is easy for 
numerical calculation. We will use wavelet D4. The basic scaling function (/)(x), basic wavelet 
tp{x) and their Fourier transform of the D4 are shown in Fig. 1. 

To simplify the notation, we consider an 1-D density field n{x) and velocity field v{x) on 
spatial range L. The result is straightforward generalized to 3-D fields in §3.3. 

For doing the DWT analysis, the space L is chopped into 2^ segments labelled by 
I = 0,1, ...2-' — 1. Each of the segments has size L/2K The index j can be a positive 
integral. It stands for scale L/2K The index / is for position, and corresponds to spatial range 
lL/2^ <x< {l + l)L/2K 

For a random field n(x), the DWT analysis is performed by the scaling functions 
(l)j^i{x) = {T /LY/'^(I){V /L - I), and wavelets ipj^x) = {T /LY/'^ip{T /L - I). The scaling functions 
play the role of window function. Generally, (pj^i{x) is a window in the segment /. They are used to 
calculate the mean field in the segment /. The wavelets tpj^i{x) essentially is similar to the function 
U of eq.(lO). They are used to extract the fluctuations of the flelds at the segment /, i.e. they are 
used to calculate the difference between the mean fields at space ranges lL/2^ < x < [l + l/2)L/2^ 
and (/ + l/2)L/2i < x < {I + l)L/2L 

The scaling functions and wavelets ipj^i^x) satisfy the orthogonal relations as 

4>j,iix)4>j,i'ix)dx = Si^i>, (12) 

ijj^i{x)'ipj>^if{x)dx = Sjj>6i^u, (13) 

(l)j,i{x)ijf^i>{x)dx = 0, if / > j. (14) 
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With these properties, a 1-D random field n{x) can be decomposed into 

n{x) = n^{x) + E E ^',l^fA^), (15) 

j'=j 1=0 

where 

™'(^) = E ^li<PjA^)- (16) 

The scahng function coefficient (SFC) e^^ and the wavelet function coefficient (WFC), e^i are 
given by 

^li = J n{x)4>j^i{x)dx, (17) 

and 

~eli = j n{x)ij,,i{x)dx, (18) 

respectively. The SFC e"^ measures the mean n{x) in the segment /, while the WFC e^^ measures 
the fluctuations of field n{x) at I on scale j. 

The first term on the r.h.s. eq.(15), n^{x), is the field n(x) smoothed on the scale j, while 
the second term contains all information on scales > j. Because of the orthogonal relation 
eq.(14), the decomposition between the scales of < j (first term) and > j (second term) in eq.(15) 
is unambiguous. The value of j in eq.(15) can be any integer, and therefore, a scale-by-scale 
decomposition becomes possible. 

Since scale r, position x, and window function W{\x' — x\/r) of eq.(3) correspond to, 
respectively, j, I, and 4>j^i{x') of the DWT analysis, the DWT counterpart of eq.(3) (for 1-D) is 

^ _ ^ v{x)n{x)(l)j^i{x)dx _ ell 
■'' J n{x)(f)j^i{x)dx e^i 

where e^; and e"; are respectively, the SFC of field v{x)n{x) and n(x), i.e. 

Sji = / v{x)n{x)(l)j^l{x)dx = Y^ WnV{Xn)4>j,l{Xn) (20) 



and 

Ng 

^1,1 = / n{x)(l)j^l{x)dx = ^ Wn<Pj,l{Xn)- (21) 



ra=l 



- 10- 

Vj^i is the mean velocity in the spatial range lL/2^ < x < {I + 1)L/2K 

Similarly, Ur corresponds to tpj^i, and therefore, the DWT counterpart of eq.(5) is 

^^ _ j v{x)n{x)'4)j^i{x)dx _ ell 
■'' J n{x)<f>j^i{x)dx e";' 

where e^; is the WFC of field v{x)n{x), i.e. 

Ng 

eji = / v{x)n{x)'ll)j^i{x)dx = ^ WnV{Xn)'^j,l{Xn)- (23) 

•^ n=l 

^Vj^i is the difference between the mean velocities of spatial ranges lL/2^ < x < {1 + l/2)L/2^ and 
{I + l/2)L/V <x<{l + 1)L/2K 

Vj^i and Avj^i are the variables of the velocity field v{x) in the DWT representation. 
These variables give a complete description of the field v{x) without loss of information. The 
orthogonality of scaling functions and wavelets insure that the decomposition does not cause false 
correlation among these variables. 

It has been pointed out that the galaxy pairwise velocity dispersion measured by conventional 
techniques, i.e. pair-counting statistics, is biased by densest regions, as the statistic actually is 
pair-weighted (Strauss, Ostriker & Cen 1998.) In the language of the velocity field decomposition, 
the conventional method is redundant for modes at the dense regions. In the DWT decomposition, 
each mode (j, /) corresponds to a cell in phase space (scale L/2^ and position I), and is described 
by one variable Vj^i or Avj^i. That is, all modes are treated equal. Each degree of freedom of 
the velocity field is represented by one variable, regardless the number of galaxies in the cell. 
Therefore, the description of Vj^i and Avj^i is free from the bias of pair-weight. 

3. The VD and PVD spectrum of a velocity field 
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3.1. The VD and PVD spectrum 

With eqs.(19) and (22), the VD and PVD are given by ensemble averages of v'^ i and Af?;, i.e. 

2v 1/2 



and 



V 12 \l/2 



<l = (Hl)'^' 



J v{x)n{x)(/)j^l{x)dx 
J n{x)(l)j^l{x)dx 



J v{x)n{x)'ipj^l{x)dx 
J n{x)(j)j^l{x)dx 



n2v 1/2 



(24) 



(25) 



If the field is statistically uniform, a^ ^ and a^ ^ are independent of /. One can define the VD and 
PVD on scale j by a'j = a^ i, and a^^ = a^^, which are referred to as the VD and PVD spectrum 
of the velocity field, respectively. 

If the "fair sample hypothesis" (Peebles 1980) holds, the ensemble average (...) can be replaced 
by a spatial average over /. Eqs.(24) and (25) become 



c,- 



-. 2^-1 



23 



1=0 



1/2 



■ . 2^-1 / V \ 2' 

-T l^l 



1/2 



and 



a- 



pv 



1 



2-'-l 



23 



J E Hi 



1=0 



1/2 



2J-1 /^v ' 



E 



/=o \^j,i, 



1/2 



(26) 



(27) 



Eqs.(26) and (27) show that the VD and PVD of a velocity field are given by the SFCs and WFCs 
of the DWT variables, respectively. For a non-gaussian random field F{x), the statistical behaviors 
of the field quantity F{x) and its spatial difference F{x) — F{x + r) generally are different. The 
statistical behaviors of the VD and PVD are different from each other. 



As a useful variant of eq.(24), we define a modified VD by 



O",- 



{[J v{x)n{x)(j)j^l{x)dxY 



1/2 



(28) 



{[Jn{x)(pj^l{x)dx]'^) 

Obviously, if the distribution n{x) is randomly uniform, a'^ will be the same as 0"^. If velocity 
field is long-range correlated, i.e. {v{x)v{x')) ~ const, and no correlation with n{x), we have also 
a^ — a'j, even n{x) is not uniform. 
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Because {[J v{x)n{x)<t>j^i{x)dxY) = (1/2-') E(ej,i) and {[J n{x)(t>j^i{x)dxY) = (l/2-')E(ej 



eq.(28) yields 



O",- 



Ei 



2^-1/1, \2 



=0 vej,i; 






1/2 



^n \2 



(29) 



Similarly, a modified PVD can be defined by 



a 



pv 



E 



2^-1 /-I, ^2 



E 



23-l|',n \2 



1/2 



(30) 



The difference between aj and aj, or a^^ and a^^ can be used as an indicator of the correlations 



between velocity and density fields (§5). 



3.2. An example: random fields 



To demonstrate the VD and PVD spectra, let's consider 1-D random fields, which are 
produced by uniformly random distribution of N particles in range L. Each particle n is assigned 
a velocity Vn drawn from a Gaussian distribution of the velocity with zero mean and variance 300 
km s~^. The two-point correlation function of number density {n{x)n{x')) = n^ = {N/ L)"^ . The 
two-point velocity correlation function of velocity field i'"{x — x') = {v{x)v{x')) would be 

(300)2 y.^2 g-2 1^ _ ^/| ^ 2./iV 

\x- x'\ > L/N . 



S,"{x — x) 



(31) 



In this case, a'j ~ a^ . Figure 1(a) displays the modified VD spectrum ct^ of the random 
samples with L = 1 and N = 10^, 10^ and 10^. Fig. 1(a) shows that a^ increases with j, and 
saturates at 300 km/s. This behavior can be deduced analytically. In fact, following eq.(28) and 
{n{x)n{x')) = const, we have 



{^D 



V.2 _ J J {n{x)n{x')){v{x)v{x'))(pj^i{x)(j)j^i{x')dxdx' 
■' J J {n{x)n{x'))(l)j^i{x)4>j^i{x')dxdx' 

! ! C{x — x')(l)j^i{x)(j)j^i{x')dxdx' 



(32) 



/ / (t)j^i{x)(t)j^i{x')dxdx' 
Considering that both (j)j^i{x) and (j)j^i{x') are window- like functions in the spatial range 
lL/2^ < x,x' < (I + l)L/2^ , the integral of the denominator of eq.(31) is proportional to area 
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L/2^ X L/2K On the other hand, the numerator is proportional to area L/2^ x L/N, the factor 
L/N is due to the condition |x — x'\ ~ L/N in eq.(31). Thus, the VD of eq.(32) is 

f 2^/2300/\/]V km^ s'^ if j < ln2 N 
a] ~ <^ (33) 

[ 300 km^ s-2 if j ~ ln2 N 

This is just what as indicated in Fig. 1(a). 

Figure 1(b) gives the modified PVD spectrum calculated by eq.(30). It shows the same as the 
VD spectrum, i.e. ct^*' = aj. This is expected. For the VD eq.(29), the SFC e" i is given by the 
addition of the velocities of galaxies within the cell {j,l), while for the PVD eq.(30), the WFC 
e^ I is given by the difference between the velocities of galaxies in Ll/2^^^ and L{1 + l/2)/2-'+^. 
However, the one point function of the velocity is Gaussian, the probability of positive and 
negative velocity is the same, and therefore, the statistical properties of the SFC and WFC are 
the same. Consequently, the VD and PVD spectra coincide for this random field. 

From Figs. 1(a) and 1(b), we can conclude that if the velocities of galaxies are uncorrelated on 
scales larger than j, the VD and PVD spectra will decrease systematically with the j decreasing 
by a law 2-''^ (1-D). This result is useful to estimate the scale of the correlation length of velocity 
field of galaxies. 



3.3. 3-D velocity field 

For 3-D velocity field v(x), x = (x^ , x'^ , x"^) , the DWT decomposition is based on the 3-D 
scaling functions and wavelets, which can be constructed by a direct product of 1-D wavelet basis 
(Fang & Thews 1998). For instance, in a 3-D volume Li x L2 x L3, we have scaling function as 

•^j.Kx) = (l)j,^l,{xi)(pj^^i^{x2)(l)j.,,Ux3), (34) 

and wavelets as 

^j,i(x) =i^j,,i,{xi)i;j2,i2i^2)i^j3,i3i^3), (35) 

where j = iJij2,J3), 1 = {h,h,h), and k = 0...2J'''~\ i = 1,2,3. The wavelet ^j,i(x) 
is non-zero mainly in a volume Li/2^'^ x L2/2^'^ x Lj,l2^'^ ^ and around the position 
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(xi = /iLi/2-'i,X2 = /2-^2/2-'^)2;3 = l^L^/2^'-'')^ i.e. they localized in both scale and physical space. 



The 3-D generalization of the DWT variables Vj^i [eq.(19)] and Aw^y [eq.(22)] are given by 



j.i 

^i j,i = 3r 



and 



An 



ij,i 






where e^-\ and ej*[ are the SFCs as 

•^ n=l 

/ ra(x)0j,i(x)dx = ^ ti;„(/)j,i(x„), 



iV„ 



«j,l 



n=l 



and e^', is the WFC as 



^jV ^ / ^« W^j.iW*^^ = X! ti;„Wj(x„)V'j,i(x„). 

n=l 



(36) 

(37) 

(38) 
(39) 

(40) 



The SFCs ej*) [eq.(39)] are assigned at regular grids U = 0...2''' and i = 1, 2, 3. That is, e]*; is an 
assignment of the number density distribution n(x) onto grids (j,l). Eq.(39) shows that the nth 
galaxy (or particle) is assigned onto grid j,l by weight (;^j i(x„) (Fang & Feng 2000). Similarly, 
SFCs e^'j is an assignment of the distribution t;j(v)n(x) onto grid j,l. Therefore, eqs.(36) and (37) 
essentially are the same as the DWT mass assignment (Fang &: Feng, 2000). We may call it the 
DWT velocity assignment. 



From eqs.(36) and (37), we have 



1 



2J-1 /,-"» \ 2' 



for 3-D VD and modified VD, and 

^^--[((A.,j,i)2)]V2 



2(ii+i2+j3) 



Tt^i& 



2J-l^gn ^2 



n 



=0 V-j,l7 



E 

1=0 V^j.i. 

1/2 



1/2 



(41) 

(42) 



O"; 



-I 2J-1 /2t'« ^ 

I V [^ 

2(ii+i2+i3) ^^ \ eV-, 

1=0 \ J,l/ 



1/2 



(43) 



-r 
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(44) 



for 3-D PVD and modified PVD. 



4. Velocity field of simulation samples 

4.1. Simulation samples 

To demonstrate tlie velocity field of cosmic matter, we produce ten realizations of the N-body 
simulation for each model of the SCDM, rCDM and ACDM. The parameters (ilo, A,r, erg) are 
taken to be (1.0,0.0,0.5,0.55) for standard CDM model (SCDM), (1.0,0.0,0.25,0.55) for a variant 
of SCDM model (rCDM) and (0.3,0.7,0.21,0.85) for a low density flat model (ACDM). We use 
modified AP^M code (Couchman, 1991) to evolve 128^ cold dark matter particles in a periodic 
cube of side length L= 256 h^^Mpc. The linear power spectrum is taken from the fitting formula 
given in Bardeen et.al. (1986). 

In our simulation, the so-called "glass" configuration is used to generate the unperturbed 
uniform distribution of particles, and the Zel'dovich approximation is then applied to set up the 
initial perturbation. For the mass assignment on the grid and the calculation of the force on a 
given particle from interpolation of the grid values, we use the triangular-shaped cloud (TSC) 
method. The starting redshift of the simulation is taken to be Zj = 15 for the SCDM model and 
Zi = 25 for ACDM and rCDM. The force softening parameter r/ in the comoving system decreases 
with time as ry oc l/a{t). Its initial value is taken to be ry = 384 h^^ kpc and the minimum value 
to be r/min = 128 h^-*^ kpc, corresponding to 15% and 5% of the grid size respectively. We use the 
"leap-frog" scheme in the single-step integration of time evolution, and take 600 total integration 
steps down to z = 0. 
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4.2. The VD and PVD spectra of the simulation samples 

Now we calculate the VD spectrum of the simulation samples by eqs.(41) and (42), and the 
PVD spectrum by eqs.(43) and (44). A problem in using eqs.(41) and (43) is the ambiguity 
caused by extremely low values of galaxy number density in voids or underdense regions in the 
galaxy distribution, i.e. n(x) ~ 0. Within these areas, the corresponding SFCs e^-\ would be very 
small, even zero, and hence, lead to very large terms of Cn/en- These terms are unacceptable, 
because in these areas the number density n(x) actually are dominated by Poisson fluctuations. 
To control these terms, we use the so-called "off-counting" algorithm, of which the details are 
given in Appendix. 

Since the cosmic velocity field is isotropic, the VD or PVD for the three components of v* and 
pv"^ i = 1,2,3 should be the same, i.e. 

a^^ = aV^ = aV^^ , af^ = ap = a^^ . (45) 

Therefore, it is enough to calculate one component among Vi and pvi. 

We consider diagonal mode, i.e. ji = J2 = js = j- In this case, the VD, a"\ and PVD, 
(T^\ depend only on one scale parameter j. The physical scale of mode j is 256/2-' h^^ Mpc, or 
wavenumber k = 2-'+^-v/37r/256. Figures 3 and 4 display, respectively, the j- or /c-dependence of 
o"^' and a^^ for models of SCDM, rCDM and ACDM. In comparison, the figures also plot the VD 
spectra for two types of random samples, which are given by assigning Ng velocities to randomly 
distributed particles (hereafter randoml), and randomizing the assignment of the Ng velocities 
with particles x„ (n = l...Ng) (hereafter random2). 

Figures 5 and 6 give, respectively, the j- or /c-dependence of a^^' and ct^^' of the simulation 
and random samples for models of SCDM, rCDM and ACDM. 

As mentioned in §3.2, if the peculiar velocities are not correlated on scales larger than j, or 
^^{x — x') ~ for |x — x'l > L/2^ , the VD and PVD spectra will decrease with the decreasing j 
according to the law of 2-' " (1-D) [eq.(33)], or 2'^^''^ (3-D). The spectra of random samples of Figs. 
3 and 5 are decreasing with j as the 2^-'' ^ law. However, the VD spectra of the simulation samples 
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are almost flat, i.e. scale-independent, in the range of A; > 0.2 h Mpc^^, and slightly decreasing for 
k < 0.2 h Mpc^^. This result shows that the peculiar velocities are spatially correlated on scales 
at least few 10 h~^ Mpc. This is probably caused by the in-fall process of forming halos on scales 
10 and 20 h^^ Mpc (Xu, Fang & Wu 2000.) 

The PVD spectrum decreases with j on scale /c ~ 1 h Mpc^"*^, but it is slower than the 2^-''^ 
law. Accordingly, the correlation of the PVD on large scales is also significant. 

Recall that in the recovery of the real space Fourier power spectrum from galaxy redshift 
surveys, to account for the redshift distortion effect, it is usually assumed that the peculiar velocity 
dispersion a^ is scale-independent. The DWT analysis show that the assumption of a constant a^ 
is reasonable on scales k > 0.2 h Mpc^^, but may have uncertainty of factor 2 or 3 on large scales. 

Although the VD and PVD of CDM models have similar scale-dependence, a difference 
between the VD and PVD can already be seen from Figs. 3-6. Figs. 3 and 4 show that for the 
CDM models, the modified VD a^^ is almost the same as a^' on all scales, while their o"^^' are 
generally larger than cjy^% especially on small scales (Figs. 5 and 6). This result indicates that 
the peculiar velocity possesses long self-correlation length, but is less correlated with local number 
density. 

This conclusion can also be seen from the behavior of randoml and random2 samples. The 
only difference between randoml and random2 is that the former has a uniform distribution n(x) 
while the latter is inhomogeneous. Figs. 3-4 shows that for randoml the modified VD is the 
same as the VD, but for random2, the modified VD is substantially lower than the VD on small 
scales. The latter is due to the clustering of density field n(x). The summation ^j X]i=o [l/^n]^ 
general is larger than l/(^ X]i=o [^n]^) ^^^ inhomogeneous distribution n(x). Thus, the result of 
cjy' ~ cjy' for CDM models means that the correlation of peculiar velocity field is not affected by 
the inhomogeneity density field. 
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4.3. The VD-to-PVD ratio (Mach number) 

The difference between the VD and PVD can be measured by the VD-to-PVD ratio as 

Mj = ^. (46) 

This ratio is essentially to measure the relative strength of the mean velocity and the fluctuations 
of velocity fleld on scale j, i.e. the ratio between the bulk velocity and velocity fluctuations. 
Therefore, it is similar to the so-called Mach number (e.g. Suto & Fujita 1990). 

For the randomized samples, the VD and PVD spectra are identical (§3.2), as both spectra 
follow the 2^-^/^ (3-D) law. Upon the deflnition of eq.(46), the Mach numbers of the random 
samples of §3.2 are equal to 1 on all scales. For simulation samples, the VD spectrum is no longer 
the same as PVD spectrum on all scales. The deviation of Mj from one is a measure of evolved 
velocity fleld. Generally, c^* is larger than o"^^' on large scales, i.e. the bulk flow is significantly 
larger than the PVD. This is due to the in-fall evolution. But a^^ might be is less than a^^^ on 
small scales. 

Figure 7 plots the Mach numbers on scales j = 1 — 7 for the three dark matter models. 
The figure shows that Mj has a peak at A; ~ 0.2 h Mpc~^ for SCDM and rCDM, and /c ~ 0.1 
for the ACDM. The peak characterizes a typical scale on which the matter undergoes the in-fall 
evolution. The clustering in the ACDM model is earlier than the other two models, and therefore, 
the ACDM's peak is on larger scale than the two others. On scales much larger than the peak, the 
mass and velocity fields are only slightly different from a Gaussian field, the VD and PVD are of 
the same order. On small scales, Mj becomes small, and finally < 1, i.e. the PVD is larger than 
the VD. This scale-dependence can clearly be seen with the ACDM model. 

It has to be pointed out that the Mach number Mj is a statistical measure of the entire 
velocity field, rather for individual object. Since the VD and PVD have different distribution 
functions (§4.5), the Mach number measured from individual objects might be significantly 
different from the Mach number of the entire velocity fleld. 
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4.4. Comparing DWT method with conventional techniques 

It would be interesting to compare the DWT measurement of VD and PVD with conventional 
techniques. Usually, the bulk velocity is given by the amplitude of the i?s-smoothed velocity field 
Vi{x) over a volume defined by a normalized window function VFr(x) of a characteristic scale R 

ViiR) = J WR{x)v^{x)dx. (47) 

We calculate the variance (T^^(i?))^'^ with gaussian window function Wr oc exp(— r^/2i?^) for 
SCDM simulation sample. The result is shown in Fig. 8. We also plot a'"'- (SCDM) in Fig. 8. 
The data a^' (SCDM) is the same as Fig. 3. Since the effective radius of the gaussian window is 
R = {3V2^y/^Rs, it corresponds to j given hy R = {3/4-iTy/^256/2^ h"! Mpc. 

Figure 8 shows that the DWT gives exactly the same variance of bulk velocity as that 
measured by gaussian smooth. 

To compare with the PDV of Fig. 5, we calculate the pairwise velocity of particles, 
Vi2{r = |r|) = |f • [v(x + r) — v(x)]|, by pair-counting method. The variance of pairwise velocity, 
(^12(^))^ ' for the SCDM sample is shown in Fig. 9. The data of o"^^' of Fig. 5 is also shown in 
Fig. 9. It is interested to see that scales less than 5 h~"^ Mpc, we have ('Ui2(^))^ — v3o"^^' and 
i^rzi^))^ — V2\/^{Vj^{R = r))^''^. This result implies that velocity distribution on these scales 
is virialized, or quasi- virialized (Xu, Fang & Wu 2000). 

On scales larger than 5 h^"*^ Mpc, (wi2(^))^ generally is larger than \/3a^^% but less than 
^/2\/^{V^ (R))^'"^ . The large value of ('Ui2(^))"'^'^ probably is because pair-counting measurement is 
not orthogonal with the Vi{R) measurement, i.e. the pair-counting techniques may be biased by 

The difference between field variable and conventional description of velocity field can easily 
be seen in linear regime of the cosmic clustering. It is well known that the linear relation between 
the density fluctuations and peculiar velocity fleld 5{x.) = —{I/HqP)^ ■ v(x) yields 

(^'(^)) = ^y^ dkP{k)WUk), (48) 
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where (y^(i?)) = J2i=i{^^ (R)) ^ P is the redshift distortion parameters, P{k) the Fourier power 
spectrum of 5(x), and Wji{k) the Fourier transform of Wii{x.). In the DWT representation, the 
(5(x) — v(x) hnear relation gives 

tt2o2 oo , 

[^^']' = (2nk ^ ^j—j—^Pinmm/2nkn2/2nkn3/2n\', (49) 

^ ' ni,n2,n3=—oo 12 3 

where [cj^]^ = X]i=i[^^']^) the Fourier power spectrum P{n) is a function of n related to the 
wavenumber hy k = 27rn/L, in? = n\ + 712 + n^, and (j) is the Fourier transform of the basic scaling 
function 0(n) (Fig. 1). 

Therefore, the only difference of the DWT measurement of VD from conventional method is 
to replace the ordinary window function W^^^k) by the DWT window \(l){ni/2^)<j){n2/2^)(l){n^/2^)\'^ . 
Fig.l shows that |(/)(n)p is similar to ordinary window. This is why {V'^{R)) ~ [o"^]^ if the windows 
have the same size. 

With the linear relation between (5(x) and v(x), the velocity dispersion within a patch of the 
velocity field sometimes is described by 

{AV\R)) = ^J^ dkP{k)[l-WUk)]. (50) 

The DWT counterpart of eq.(50) is given by the pairwise velocity. It is 
Hif3^ ^ 1 



CF' 



E . . ^ P{nmm/2^)4>{n2/2^Mm/2^)\^ (51) 



(27r)2L ^^ ni + ni + n 

^ ' ni,n2,n3=— 00 i ^ 
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where [a^^]^ = X]i=i['^?"']^i ^^"^ V' is the Fourier transform of basic wavelet i^). 

It is clearly to see the difference between eqs.(50) and (51). The r.h.s. of eq.(50) contains all 
contribution of P(k) on scale kR > 1. On the other hand, Fig. 1 shows that |^(n)p is localized at 
n ~ ±1, and therefore, the r.h.s. of eq.(51) contains only few terms around n ~ v32-'. This shows 
that eq.(51) is based on unambiguous scale-decomposition. 

As for the conventional pair-counting measurement of PVD, one cannot write down a linear 
relation between (t'i2('^))^ and P{k) as eqs.(50) or (51). Although r is used as a scale-indicator, 
the pair-counting method rely on a scale-decomposition of the distance r of pairs, but not a scale 
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decomposition of the field. Therefore, the statistics of (V^i2(^))^ ^ire different from [a^*']^ based 
on ensemble of the field variables. 

Eqs.(49) and (51) shows also the difference between the VD, Cy, and PVD, a^^ . The 
VD contains the contributions of P(n) on all scales n < 2K That is, even for a small scale 
L/2-', a^A is mainly determined by density perturbation on scales larger tha L/2K The PVD, 
cj^*', however is determined by P{n) on the scales L/2K From the DWT analysis, we have 
|V'(n.)p = \(f){n/2)\^ — |(/>(n)p. Therefore, the PVD is determined by the power of density 
perturbations in the scale range from L/2^ to L/2^^^. 

Generally, density perturbations on large scale is linear, while non-linear on small scales. 
Therefore, eq.(49) would be a good approximation for the VD, even on small scales, while eq.(51) 
would be a poor approximation for the PVD on scales for which the non- linear clustering is onset. 
This can be seen from Figs. 8 and 9. For the VD (Fig. 8), the theory [eq.(49)] and numerical 
simulation gives basically coincident results. However, for the PVD (Fig. 9), the result given by 
eq.(51) is significantly lower than numerical simulation. Therefore, one may conclude that the VD 
and PVD describes, respectively, the linear (large scales) and non-linear (small scales) behavior of 
the cosmic mass field. 



4.5. One-point distribution of Vi and Awj 

The PDF of the galaxy pair velocity is usually modeled as an exponential, 
f{Av) oc e~^ Ad/o-p"^ This model is supported by best fitting of early galaxy redshift 
surveys, such as 14.5 m^ CfA (Davis & Peebles 1993) and 1.2 Jy IRAS (Fisher et al. 1994.) 
Some non-linear clustering models also yield exponential distribution of pair velocity (Sheth 1996, 
Diaferio & Geller, 1996). However, these results were measured by conventional techniques, it 
cannot directly be applied to variables Vi jj or Avi jj. We calculate one-point distributions of field 
variables Vi and Awj. 

For a given scale j, the 2'-'^"'"-^2+i3) values of vi j\ (or Cn/ej,!) form an ensembles of Vi j\. Thus, 
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the one-point distributions f{vi) on scale j can be obtained directly from the distribution of Vi jj. 

Figure 10 gives the one-point distribution f{vi) on scales j = 2 — 7 of the SCDM. Models 
of rCDM and ACDM give the similar results. Figure 10 shows that the distributions are 
scale-dependent. On large scales j = 2,3 the distribution is close to a Gaussian. On small scales 
j = 6, 7 the distribution in the middle range {\vi\ < 450 km/s) can still be fitted with a Gaussian 
function. In the range 400 < Vi < 1200 km/s, f{vi), however, follows a straight line which implies 
that the distribution is exponential. With increasing Vi > 1200 km/s, f{vi) decays slower than a 
straight line. It indicates that the tail of the one-point distribution f{vi) extends further than an 
exponential distribution. Therefore, only in the middle velocity range, we can fit the distribution 
f{vi) by an exponential. 

From Fig. 8 and eqs.(48) and (49), one can expect that the PDF of VD measured by a 
conventional Gaussian window should be about the same as Fig. 10. As has been mentioned, the 
scaling function of the DWT analysis is just a window function. It will give the similar statistical 
result as a Gaussian window if the orthogonality and completeness are not the key of the relevant 
statistics. 

Figure 11 gives the one-point distribution f(Avi) on scales j = 2 — 7 of the SCDM. Models 
of rCDM and ACDM also give the similar results. The one-point distributions f{Avi) are also 
scale-dependent, but it is significantly different from f{vi). On large scale j = 2,3 the distribution 
f{Avi) is approximately exponential, i.e. ln/(Awj) vs. Avi can roughly be approximated by 
a straight line. However, on all scales j > 5, the distribution f(Avi) is nothing but typically 
lognormal, i.e. ln[f (Avi) / f (2'^ Avi)] = a\n[f (Avi) / f {2Avi)], where a is a real number. Lognormal 
random field is often employed to model the non- linear clustering (Cole &: Jones, 1991). The 
mass field traced by the Lya forests is also lognormal (Bi & Davidsen, 1997; Feng & Fang 2000). 
Moreover, lognormal random field is typically intermittent (e.g. Zel'dovich, Ruzmaikin, & Sokoloff 
1990). Based on these considerations, it is not unusual to find that the one-point distribution of 
Avi is lognormal. 

We also calculate the one-point distribution of pairwise velocity measured by conventional pair 
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counting. The result is ploted in Fig. 12, which are the same as all measurements by pair counting 
method, i.e. the distributions basically are exponential. Although we have {vi2{r))^''^ — \/3<7?"' 
on small scales, the PDF of Vi2{r) is very different from At;j in the same scale range. 

5. Correlations between velocity and density fields 

5.1. Local velocity-density correlation 

In redshift distortion theory, VD, a"" , and its distribution generally are assumed to be 
independent of the density field. Observational evidences, however, seem to support the existence 
of correlation between the velocity and local number density of galaxies. For instance, the pairwise 
velocity dispersion measured by the conventional method is found to be sensitive to the presence 
of dense objects, like rich clusters (e.g. Mo, Jing & Borner 1993). However, the conventional 
method of measuring the peculiar or pairwise velocity is biased to dense objects, due to the 
number-counting and pair-counting-weighted statistic. Thus, the f (x) — n(x) correlation might be 
contaminated by the density-density n(x) — n(x) correlation. With the field variables vi j\ and 
Avj j^i, one can distinguish among the correlations of w(x) — n(x) and n(x) — n(x). Thus, one can 
detect the correlation between the velocity and local density without bias. 



Let us consider the number-counting statistics of the mean peculiar velocity. It is given by 



.. ... = ^-V"""-^""^ (52) 



Using the "partition of unity" of wavelet (Daubechies 1992), i.e. X]i=o </*] i(^) — 
{LiL2L3)/2^''-^^^^^'-^ , one can rewrite eq.(52) as 

l^i trad — ^i-i N / x ~ 2J-1 ' V"-*"^/ 

El=0 E„=l ^n</'j,l(Xn) El=0 ej"l 



Therefore, the existence of the local v-n correlation can be tested by 



j trad Y^2J-1 1^1)^12 Y^2J-1 i^n 12' ^ ' 



Y^z;j-i \ Vi |2 sr^-z-i-i \ r, 



j,ii 
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where we still consider the diagonal modes j'l = j2 = Js = j. If there is no correlation between the 
local velocity variable e^'j and density variable ej^j, we have C* ^^'"^ = 1- We display the result 
in Fig. 13, which shows that for the CDM model C- ^^^^ ^ 1 on A; > 1 h Mpc^ . Nevertheless, 
we should not immediately conclude the existence of strong correlation between f (x) and n(x), 
because Fig. 13 also shows a strong deviation from C- ^^^^ = 1 for random2, i.e. the samples 
without velocity-density correlation. That is, the C- ^^^^ ^ 1 given by random2 is completely 
from density-density correlation. 

To identify the v-n correlation, it is essentially to refer to both randoml and random2. The 
values of C- ^^^^ for the CDM mode deviate not only from randoml, but also random2. Thus, we 
may conclude that the peculiar velocity given by the number-counting statistics is correlated with 
density field, but not very strong. 

For the DWT velocity variables f j j,i [eqs.(19) and (36)], the mean peculiar velocity is 
(1/2J) X]i=o ^Ti/^n- Thus, the local v — n correlation can be measured by 






That is, the ratio between a"^ and &"a is a measure of the VD - n correlation. The CDM spectra of 
0"^ and fj^ shown in Figs. 3 and 4 are very similar. Therefore, in term of the DWT variables, the 
VD-n correlation is very weak (see Fig. 13). 

The traditional estimation of the velocity dispersion within the patches (j, 1) is given by 

X^^Vi trad\ = ^nj 1 • (56) 

^1=0 *"],! 

It actually is the mean value of the absolute pairwise velocity of models (j, 1). The local correlation 
between this pairwise velocity and density field can be estimated by 



Y^2J-l |^«i|2 Y^2J 

Z^1=0 l^j.ll 2.^1=0 l'=j,ll 



j trad Y-2J-1 |~«i|2Y-2J-l I „ |2' ^ ^ 



The result of C ■ ^^^^ is given in Fig. 14. Although for the CDM model, the correlation C- ^^^^ is 
much larger than 1 on small scales, but it is always less than the C- ^^^^ of random2. Therefore, 
it seems to be not an evidence of the correlation between pairwise velocity and local density. 



25 



Using the field-variable-defined pairwise velocity Atij j i, the correlation should be can be 
measured by 



2^1=0 l^^»j,il 2^1=0 l^j,il ^ ■? 



Fig 14 also displays this ratio. For the CDM models, C- ' deviates from C- =1 

(randoml), but it is similar to that of random2. Therefore, no strong correlation between the 
pairwise velocity and local density can be identified. 



5.2. Correlation between PVD and local density fluctuations 

To account for the redshift distortion effect in measuring the DWT power spectrum, it would 
be necessary to investigate the correlation between the density fluctuations (5(x) and velocity 
field v(x). The DWT power spectrum is determined by the density fluctuations on mode (j,l). 
The bulk velocity of mode (j,l) causes only a position shift of mode (j,l) in redshift space with 
respect to real space, but do not change the corresponding power of the density fluctuations of 
mode (j, 1). On the other hand, the pairwise velocity of mode (j, 1) will change the power of the 
density fluctuations of mode (j, 1). Therefore, the redshift distortion of the DWT power spectrum 
is mainly caused by velocity field variables Avj i, not Vj i. We will only study the correlation 
between Avj i and density fluctuation 6(x). 

A simplest statistical measure of the Av; i-6 correlation is the second order statistic defined as 



^2J — 1 ~v^ ~n 



.{pv,S) _ ^ 2.^1=0 ^j,l^j,l 

E2J — 1 |~^i| Y->2J — 1 \~n I 
1=0 l^j.ll Z^1=0 l^j.ll 



3 v^2J-l \~vi\ v^2J-l 



(jyv^^o) -^viH be equal to zero, if there is no correlation between velocity fluctuation WFC ("A and 
density fluctuation WFC l^y 

Fig. 15 shows C- ^' in a CDM model, all results are consistent with (7- ' = within 
the error bars. Thereby, in terms of second order statistics, it is likely reasonable assumption 
that there is no Avj^i-(5 correlation. However, C- ' = doesn't mean that Avj^ and (5(x) are 
uncorrelated in general. Since the one-point distribution of Avj i are invariant with respect to 
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Avi{x) ^ —Avi{x), in this case, statistics with hnear term of Avj(x), hke (Afj(x)(5(x)) wih 

be zero. Therefore, we should study the correlation between |At;j(x)p and 5^(x), which can be 

calculated by 

Mpv\s^) ^ ^' 2^1=0 [gj,igj,iJ .„„. 

j ^2J-l |~t,M2 Y-2J-1 |~n |2' ^ ^ 

The results for the dark matter model are also plotted in Fig. 15. It shows that C- ' of the 
CDM models is higher than that of randoml and random2. Therefore, Av and (5(x) are locally 
correlated on scales k > 0.2 h Mpc~^. 



6. Conclusion 

It is known that in the DWT representation a continuous field, like cosmic density field n(x), 
can be described by discrete variables, given by the DWT scale-space decomposition of n(x) into 
the DWT mode (j, 1) (Fang & Feng 2000). Since the modes (j, I) are complete and orthogonal, the 
DWT decomposition is information-lossless and not redundant. It insures that one can compare 
statistical results of different samples. 

Taking the advantage of the DWT analysis, we show that the cosmic velocity field can also be 
properly described by the discrete variables Vi j j, which are given by an assignment of the number 
density and velocity of galaxies (or particles) into the DWT modes (j,l). In this scheme, the 
peculiar velocity and pairwise velocity of galaxies or particles are given by field variables, which 
correspond to the coefficients of the scaling functions and wavelets of the DWT decomposition. 
As a consequence, the VD and PVD are no longer measured by number-counting or pair-counting 
statistic, but with the ensemble of the field variables. 

Using simulation samples, we show that peculiar velocity field is significantly different from 
randomized field. The peculiar velocities show correlation on scales of few tens h~^ Mpc. The 
pairwise velocity (or relative velocity) has similar correlation. On small scales, we also found 
significant correlations between pairwise velocity and density fiuctuations. This is especially 
valuable in treating the redshift distortion of the DWT power spectrum. 
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Although the VD and PVD look similar to each other in some aspects, they actually are very 
different both statistically and physically. The VD-to-PVD ratio shows the difference between 
the scale-dependencies of the VD and PVD. More prominent difference between the VD and 
PVD is shown by one-point distribution. The one-point distribution of the VD is approximately 
exponential, while the PVD is lognormal. This difference is typical of intermittent field (a 
similar difference between the one-point distributions of field variable A and the difference AA in 
turbulence is given by Shraiman & Sigglia 2000). Therefore, one can conclude that the cosmic 
velocity field is intermittent. 

LLF and YQC acknowledges support from the National Science Foundation of China (NSFC) 
and National Key Basic Research Science Foundation. 

A. Poisson noise 

The observed galaxy number density distribution n(x) [eq.(2)] is believed to be a realization 
of a Poisson point process with an intensity n(x) = n^(x)[l -|- (5(x)], where n^(x) is selection 
function, and (5(x) = [n(x)/n^(x)] — 1 is the density contrast fluctuations of the underlying matter 
field. Thus, the correlation function of 5{x) is 

(5(x)5(x)) = -1 + /^^M!^\ - Soi^ - xO^TT. (Al) 

The power spectrum of (5(x) is then given by (Fang & Feng 2000, Yang et al. 2001) 

Pj = /?-Arj. (A2) 

where 



^j " 20-i+i2+i3) 2^ 2^ 2^ T^02 (^3) 

/l=0 «2=0 /3=0 ^ J-1^ 



and 



2J1-1 232-1 2W-1 -. 

^j " 2(ii+i2+i3) 2^ 2^ 2^ ^- (^^) 

/l=0 «2=0 Z3=0 J>1 



28 



where the WFC e?j is given by 

el, = lni^)i^i^,i^)dx (A5) 

and n?j is the mean number density of the selection function n^(x) in the mode (j, 1). 

The term I? is the mean power of j modes measured from the observed reahzation n^(x), and 
the term A'j is the power on j modes due to the Poisson noise. Therefore, the Poisson error for the 
galaxy distribution n(x) can be described by an error distribution as 

2"'i-12"'2-12"'3-l 

E E E (A6) 

ll=0 1-2=0 l3=0 
1/2 

(5^(X1 - /iLi/2-^i)5^(x2 - /2L2/2^2)<5^(X3 - /3^3/2-^3). 



L1L2L3 _„ 

n 



2J1+J2+J3 J.i 

where (Ji, J2, J3) correspond to the smallest scale of sample considered. In the case of n^(x 
const and Li = L2 = L3 = L, Ji = J2 = J3 = J, Eq.(A6) becomes 



2^-1 ,^ 



W= E\ ^'^''(^-1^2^) (A7) 



1=0 



It simply means that the Poisson noise of the galaxy number in a cell (J, 1) is J Ng/2^-^ . 

B. OfF-counting algorithm 

The typical terms in eqs.(41) and (43) are Cn/en and en/en, where the denominator is given 
by eq.(39) 

^jV = / ^W'^j.lW^X = E ^i'n0j,l(Xn)- (Bl) 

n=l 

Therefore, for modes j, 1 corresponding to position of the voids of galaxies, i.e. n(x) is very small 
or zero, the SFC ej^j would be very small, even zero. These modes will lead to very large term 
of e^'j/ej*). Actually, these terms have very large error, because the number density n(x) in these 
modes generally is less than the error distribution at the same spatial range rf^'^'^{ii). 

To avoids the contamination of these large error term, we use the following algorithm. 

1. Perform DWT for the number density distributions of n(x), and rf'^'^^'^. 
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2. The SFC of n(x) generally is larger than the SFCs of n^'"''"'' (x) , and therefore e^^ > e'^J^"'^ . 
However, for some modes (j,l), we have 

^n ^ ^error /"DO^ 

^j,i < ej,i (B2) 

These (j,l) are noise dominated modes, and no observed data are available. 

3. We eliminate all modes affected by the noise dominated modes. This is, if (j,l) satisfies 
eq.(B2), we off-count the following modes 

Ji > n (B3) 

I'. = 12^H^....{1 + 1)2^'-^' -1. 

This is called off-counting algorithm. It is effectively to eliminate the effects of e^^ ~ 0, and 
give reasonable estimation of VD [eq.(41)] and PVD [eq.43)] on the same scale range as the power 
spectrum. 
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of Daubechies 4 wavelet. 
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Fig. 2. — The VD spectrum, ct^ [panel (a)] and the PVD spectrum, a^^ [panel (b)] of 1-D random 
samples, which are uniform random distribution of N particles in spatial range L = 1. Each 
particles assigned with a velocity Vn drawn from a Gaussian distribution with variance 300 km ■sT^ . 
i is for scale 1/2-'. The error bars are l-ir variance from 1000 realizations 
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Fig. 3. — The VD spectrum of the diagonal mode, a^' of simulation samples for models SCDM 
(upper panel), rCDM (central panel) and ACDM (lower panel). The randoml samples are given by 
assigning N velocities to randomly distributed particles, and random2 is produced by randomizing 
the assiernment of A'^ velocities to N mass narticles. The error bars are l-a variance from 10 
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Fig. 4. — The modified VD spectrum, a^* of diagonal mode (ji = J2 = js = j) for SCDM simulation 

samples (upper panel), rCDM (central panel) and ACDM (lower panel). The random samples are 
given by the randomizations as Fig. 3. The error bars are 1-a variance from 10 realizations for each 
model. 
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Fig. 5. — The PVD spectrum, cTj, of diagonal mode (j'l = J2 = Js = j) for SCDM simulation 
samples (upper panel), rCDM (central panel) and ACDM (lower panel). The random samples are 
given by the randomizations as Fig. 3. The error bars are 1-a variance from 10 realizations for 
each model. 
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Fig. 6. — The modified PVD spectrum, cjy^% of diagonal mode (ji = J2 = js = j) for SCDM 
simulation samples (upper panel), rCDM (central panel) and ACDM (lower panel). The random 
samples are given by the randomizations as Fig. 3. The error bars are 1-a variance from 10 
realizations for each model. 
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Fig. 7. — The VD-to-PVD ratio or Mach number Mj, of simulation samples for models SCDM 

(upper panel), rCDM (central panel) and ACDM (lower panel). The random samples are the same 
as Fig. 3. The error bars are 1-a variance from 10 realizations for each model. 
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Fig. 10. — The one-point distribution of Vi on scales j = 2 — 7 for simulation samples of the SCDM 
model. 
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Fig. 12. — The one-point distribution of ^;l2 measured by the traditional pair counting method for 
simulation samples of the SCDM model, r is the distance of the pairs in the unit of h^^ Mpc. 
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Fig. 13. — The fourth order correlation between the traditional estimated VD and local number 
density of particles [panel (a)], and the fourth order correlation between the VD and local number 
density of particles [panel (b)] for model SCDM. The random samples are the same as Fig. 3. The 
error bars are l-rr variance from 10 realizations. 
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Fig. 14. — The fourth order correlations CYtrad ' [^<l-(57)] (upper panel), and C ' [eq.(58)] 
(lower panel) for model SCDM. The random samples are the same as Fig. 3. The error bars are 
1-0" variance from 10 realizations. 
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Fig. 15. — The second (upper panel) and fourth (lower panel) order correlations between the PVD 
and local density fluctuation of particles for model SCDM. The random samples are the same as 
Fig. 3. The error bars are 1-a variance from 10 realizations. 



